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Abstract Synthetic aperture radar (SAR) has a pivotal role 
as a remote imaging method. Obtained by means of coherent 
illumination, SAR images are contaminated with speckle 
noise. The statistical modeling of such contamination is well 
described according with the multiplicative model and its 
implied distribution. The understanding of SAR imagery 
and scene element identification is an important objective 
in the field. In particular, reliable image contrast tools are 
sought. Aiming the proposition of new tools for evaluating 
SAR image contrast, we investigated new methods based on 
stochastic divergence. We propose several divergence mea- 
sures specifically tailored for distributed data. We also 
introduce a nonparametric approach based on the Kolmo- 
gorov-Smirnov distance for data. We devised and as- 
sessed tests based on such measures, and their performances 
were quantified according to their test sizes and powers. Us- 
ing Monte Carlo simulation, we present a robustness anal- 
ysis of test statistics and of maximum likelihood estima- 
tors for several degrees of innovative contamination. It was 
identified that the proposed tests based on triangular and 
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arithmetic-geometric measures outperformed the Kolmogo- 
rov-Smirnov methodology. 

Keywords robust statistics • information theory • nonpara- 
metric methods • parametric inference 



1 Introduction 

Synthetic aperture radar (SAR), ultrasound-B, sonar, and la- 
ser images are obtained with coherent illumination and, con- 
sequently, their appearance is affected by a signal-depen- 
dent granular noise called speckle |41|. Due to deviations 
from classical properties of additivity and Gaussian distribu- 
tion, speckled images require tailored processing and analy- 
sis techniques which could conveniently rely on the statisti- 
cal properties of the data. 

The ability to discriminate SAR image regions is an im- 
portant tool for identifying and classifying distinct targets 
in the scene. Such capability has also been employed as a 
validation criterion for novel methodologies in SAR image 
change detection |26|. These facts have motivated the pro- 
posal of contrast measures based on the statistical properties 
of image data II 41127 l[39ll . Therefore, an accurate modeling 
of speckled data is a major step for SAR image analysis ifTSl . 
Indeed, the probability distribution that describes the image 
data is a crucial information for any subsequent statistical 
analysis. 

Among the proposals for SAR image modeling, the mul- 
tiplicative model 1 13] has proven to be an important frame- 
work for the statistical processing and analysis of SAR 
data O. Assuming that the hypotheses suggested by the 
multiplicative model are valid, the distribution arises as 
an efficient and effective probabilistic model for speckled 

Hypothesis test methods have been considered as a 
venue to quantifying the image contrast between different 
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regions in SAR imagery. Edge detection ||9|[2Q||, classifica- 
tion 12411331 , target identification |23 |, change detection | 38 | 
and oil spill identification | 34 | rely on parametric and non- 
parametric hypothesis tests. Although some nonparametric 
methods are computationally convenient, they seldom take 
explicitly the roughness information into account \9\. As a 
consequence, such simple approaches are unable to provide 
any information about that important parameter. Addition- 
ally, these methods cannot offer any known statistical prop- 
erty that could be utilized for the design of hypothesis test- 
ing procedures |20|. Thus, the proposal of hypothesis testing 
methods based on roughness dependent contrast measures is 
a much sought tool for speckled data analysis. 

In recent years, the interest in adapting information- 
theoretic tools to image processing has increased notably. 
In 1 39 1, the Bhattacharyya and (symmetrized) Kullback- 
Leibler distances were applied as scalar contrast measures 
for polarimetric and interferometric SAR imagery. The 
works |[T7ll40 l proposed several parametric methods based 
on (/z, (^) -divergence measures for speckled data. 

Often regarded as a classical benchmark tool, the Kol- 
mogorov-Smirnov distance has been used for segmentation 
and labeling of SAR images |32|. This nonparametric test 
has found applications in polarimetric speckled data anal- 
ysis |24|. Additionally, the Kolmogorov-Smirnov distance 
was also applied as a goodness-of-fit measure for SAR data 
modeling O. 

Whereas on the one hand a careful statistical modeling 
is fundamental, on the other hand nature is more complex 
than most models. Even sophisticated statistical models and 
tools may have their effectiveness diminished when data de- 
viate from the assumed underlying hypotheses, even mildly. 
Robust statistics is a way of tackling this problem, and re- 
sistant techniques have been used with success in the SAR 
literature dElIIIllaa. 

This paper presents three contributions regarding the 
analysis of SAR imagery with stochastic measures of con- 
trast. Firstly, it quantifies the ability of selected paramet- 
ric methods based on divergences measures by means of 
the size and the power of associated hypothesis tests using 
Monte Carlo. Such quantities were compared to the ones re- 
lated to the nonparametric Kolmogorov-Smirnov test. It also 
presents an expression for the Kolmogorov-Smirnov dis- 
tance between distributions. Secondly, it employs a con- 
tamination model for representing 'innovative outliers' |[T2l 
and assesses the behavior of maximum likelihood (ML) es- 
timators under the model in the presence of contami- 
nation. Finally, it submits hypothesis tests to different de- 
grees of contamination, aiming to investigate the robustness 
of their sizes. Additionally, it presents an application to ac- 
tual SAR data. 

The remainder of the paper is organized as follows. Sec- 
tion discusses the multiplicative model and the dis- 



tribution. The hypothesis tests under consideration, namely 
those derived from (/z, 0) -divergences and the Kolmogorov- 
Smirnov distance, are derived in Section [3] Section |4] 
presents Monte Carlo experiments in order to assess ML es- 
timation and hypothesis tests under situations which include 
"pure" and contaminated data. Actual SAR data is analyzed 
in Section|5] Section [6] concludes the paper. 

2 The multiplicative model and the distribution 

Meaningful statistical processing of SAR imagery is always 
associated to a statistical model. 

In 12111 . Gao provided a comprehensive exposition of 
existings models. In particular, three major categories of 
models were identified: (i) the multiplicative model (131 
[161, (ii) models based on the generalized central limit the- 
orem 1 29 1, and (iii) the empirical distribution model IS. 
Other approaches include the joint distribution model lU, 
the mixed Gaussian model |11|, mixtures of Gamma 
laws 1 49 1, and the correlation based model |30|. 

The multiplicative model has a phenomenological na- 
ture which is closely tied to the physics of the image forma- 
tion Bn . Zhang BtII showed evidence that the distributions 
stemming from the product between independent random 
variables which describe the speckle noise and the backscat- 
ter outperform other models. In (21] p. 776], Gao includes 
the following note about the multiplicative model: 

The most attractive achievement among them is the 
statistical modeling on extremely heterogeneous re- 
gion of SAR images proposed by Frery et al. (T6ll 
who [. . . ] has introduced the original idea that for the 
purpose of statistical modeling, SAR images can be 
divided into homogeneous regions, heterogeneous 
regions and extremely heterogeneous regions, ac- 
cording to their contents. 

Such assertion, related to distribution, is due to the fol- 
lowing facts: 

1 . it has as many parameters as the K law, i.e., scale, rough- 
ness and number of looks, and they can also be inter- 
preted, leading to scene understanding; 

2. the roughness is expressive, in the sense that it is able to 
describe more targets than the K law (36l ; 

3. its density does not depend on special functions, as is the 
case of the K law 1 16 |; 

4. its cumulative distribution function can be obtained 
using its relationship with the Snedekor's F distribu- 
tion |37|. 

We, therefore, admit the law as the model for the ob- 
served data, regardless the target and, without loss of gener- 
ality, in its intensity format, i.e., quadratic detection. 
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Consider positive and independent random variables X 
and Y as models for the terrain backscatter and the speckle 
noise, respectively. The model for the return isZ = X Y. 

The generalized inverse Gaussian distribution is pro- 
posed as a general model for the intensity backscatter. How- 
ever, a tractable particular case is the reciprocal gamma 
law |[T3l , whose density function is given by 



fx{x;a,y) 



r{-a) 



x«-iexp(-^), 



-a,7,x>0. (1) 



In intensity format, the random variable Y obeys the 
gamma distribution with unitary mean [13 1, which has the 
following density 



fY{y;L) 



/-iexp(-L};), y>0,L>l, 



r{L) 



(2) 



where L is the number of looks. Throughout this paper, the 
number of looks is assumed to be known and constant over 
the whole image. A detailed account of the until recently 
largely unexplored issue of estimating L is provided in (31 . 

Considering independent random variables with the den- 
sities presented in Equations ([T]) and ([2]), we obtain that the 
density of Z is expressed by 



/z(z;a,7,L) : 



■ a 



7«r(-a)r(L) 
r{-a)r{L)\rJ (i. 



(3) 



where z > 0, a < is the roughness, 7 > is the scale, and 
L > 1. We indicate this situation as Z ~ ^^(a, 7,L). 

As shown in I36II37I , this distribution can be used as an 
universal model for speckled data. The rth moment of Z is 
expressed by 



r{-a-r)r{L + r) 



(4) 



r(-a) r{L) ' 

if — r > a. Otherwise, it is infinite. 

Reparametrizing Equation ([3j with M = — a > and 7 = 
M/i, we obtain the density of the Fisher model [45 J whose 
density is given by 



/z(z;M,Ai,L) : 



r(L + M) / L 



r{M)r{L)\Mii) (i + ^z) 

Using some Mellin transform properties for Fisher dis- 
tribution, Galland et ah |19| derived the so-called log- 
cummulants w/ = E(Z^ logZ): 

wo=log^ + /(L)-/(M), 



W2 ■ 



M-l 



log^ + /(L+l)-/(M-l) 



m2(L+1) 
"(M-1)(M-2)L 



log 



•V/^^(L + 2)- VA^(M-2) 



where is the digamma function. As the Fisher distribu- 
tion is a reparametrization of the distribution, the results 
in 1 19,45] apply to the latter. 

Fig. [T] shows sixteen 256 x 256 pixel regions of 
simulated speckled data according to the distribu- 
tion for L = {1,8}, a e {-15,-8,-5,-3}, and E(Z) = 
/i G {1,10}. The homogeneity depends on a, the rough- 
ness parameter, wh ile { l controls the mean bright- 
ness. Fig. 2(a) and 2(b) show the density of dis- 
tributions with varying (a,L) G {(-1.5,20), (-1.5, 8), 
(-1.5,3), (-15,20), (-15,8), (-15,3)}, and ^ G {1,2,5} 
for /i = 1 and (a,L) = (—2.5,3), respectively. Increasing L 
or /i flattens the densities, while small values of a lead to 
images of low variability. 




Fig. 1 Synthetic imagery 
{(1,1), (1,8), (10,1), (10,8)} (left to 
{-15, -8, -5, -3} (top to bottom). 
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Fig. 2 densities. 
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Parameters a and 7 can be estimated according to sev- 
eral methods, including bias-reduced procedures |[6l l43ll46ll 
and robust techniques (such as L- and M-estimators adapted 
to asymmetric distributions) dO. Small sample methods 
are particularly relevant in image processing and analysis, 
for instance, in the determination of the equivalent number 
of looks, in adaptive filters and for training classifiers. In 
such cases, samples are usually very small when compared 
to the full image 1 17], ranging from a single observation to, 
typically, less than a hundred. 

ML estimation was adopted due to its optimal asymp- 
totic properties |42|. Let Zi, . . . ,Z„ be ^ iid ^^(a, 7,L) dis- 
tributed random variables. The ML estimator for (a, 7), for 
a given L, namely (S,7), is any solution of the following 
system of non-linear equations f40l : 

VA^-a) - i/aO(L- a) -log7+ ^Itilog(7+^Z0 = 0, 

-| + %^i:Li(r+^zo-i=o. 



In general, the above system of equations does not possess 
a closed form solution and numerical optimization methods 
are required. We used the Broyden-Fletcher-Goldfarb-Sha- 
nno (BFGS) method, which is regarded as an accurate tech- 
nique Q. Next section presents several contrast measures 
for distributed images based on the above estimators. 

Fig. |3(a) shows an image over San Francisco, CA, ob- 
tained by the AIRSAR sensor at band L with four (nomi- 
nal) looks. Urban and forest regions were selected, and Ta- 
ble[T]presents the sample sizes and ML estimates of a and 7. 
Moreover, this table shows the sum of squared errors (SSE) 
between the histogram and the fitted density /„, which is 
given by 1481 



SSE: 



# pixels / ^ 



■fnf 



# pixels 



Two densities are considered: the basic gamma model (the 
scaled version of the density presented in Equation (|2])) and 
the distribution. 



Table 1 ML estimates and SSE 











SSE 




Regions 


Z 


a 


7 (xlO-2) 


^0 


# pixels 


Urban 


0.374 


-1.370 


19.550 


0.062 0.318 


2277 


Forest 


0.260 


-1.275 


7.461 


0.108 0.389 


1813 



The values of SSE provide evidence that the distri- 
bution outperforms the gamma distribution for the consid- 
ered regions. The strongest difference occurs in urban re- 
gions, i.e., in situation with intense roughness, as presented 
in Figs. |3(b)|3(c)| As discussed in (361 , among other ref- 
erences, the law includes the gamma distribution as a 
particular case and, therefore, it is able to describe homoge- 
neous targets as pasture, crops, and sea. 




(a) SAR image with selected regions 




0.5 1.0 

Gray level 



(b) Urban region 



0.2 0.4 0.6 

Gray level 

(c) Forest region 



Fig. 3 Real SAR image with selected regions and empirical and fitted 
densities for differents urban and forest regions. 



3 Contrast measures for the distribution 

An image can be interpreted as a set of regions described by 
possibly different probability laws, and contrast analysis ad- 
dresses the problem of quantifying how distinguishable two 
image regions are. Many statistical approaches have been 
developed to reach this goal, leading to the difficult problem 
of specifying accurate, expressive, and tractable metrics. In 
general terms, these methods can be either parametric or 
nonparametric. 

In a previous work, we proposed a class of parame- 
tric measures based on stochastic distances. Advancing that 
study, among the existing nonparametric methods, we now 
include an approach based on the Kolmogorov-Smirnov 
test. This method is reported to possess good discrimina- 
tory properties (311 and has found applications in several 
instances 12411321 . In the following, we briefly describe the 
parametric methods. Moreover, we adapt the Kolmogorov- 
Smirnov test to distributed data. 
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We assume that Z' and are distributed random 
variables with cumulative distribution functions F^i and Fz" , 
respectively. These distributions share the support and 
are defined over the same probability space, equipped with 
densities fz'{z\Q\) and fz"{z\Qi), respectively, where Q\ 
and 02 are parameter vectors. 



Section 3.1 presents contrast measures derived from 
stochastic distances, while Section |T2l discusses the use of 
the Kolmogorov-Smirnov distance for the model. 



3.1 Contrast measures based on divergence 

Information theoretical tools known as divergence measu- 
res offer entropy based methods for contrasting stochastic 
distributions. Divergence measures were submitted to a sys- 
tematic and comprehensive treatment and, as a result, the 
class of -divergences was proposed |42|. 

The (/z, 0) -divergence between f^i and fz" with common 
support / is defined by 

where : (0,oo) [0,oo) is a convex function, h\ (0,oo) 
[0,oo) is a strictly increasing function with /z(0) = 0, and 
indeterminate forms are assigned value zero. Table |2] shows 
the h and functions considered in this work. 

These divergences are often asymmetric, and in order 
to solve this issue and to turn them into distances, a new 
measure can expressed as 



4(z',z'0 = 



D\(Z^,Z'')^D\(Z^\Z^) 



(5) 



We will be interested in comparing samples from pos- 
sibly the same distribution ^(0), so instead of indexing 
Equation ([Sj by the random variables, we will do it by the 
parameters 0i and 02- Furthermore, denote 0i and 02 the 
maximum likelihood estimators of 0i and 02 based on in- 
dependent samples of sizes m and n, respectively. 

Lemma 1 Let the regularity conditions proposed in 
p. 380] hold. If -f- > X e (0,1) and 61 = 02, then 

m-\-n yy. fJ^OO 



2mn 



1 



m + ^/z'(0)f'(l) 



^ 2 



where 



denotes convergence in distribution. 



Based on Lemma [T] the following statistical hypothesis 
test for the null hypothesis 0i = 02 can be derived: 



5'1(01,02 



2mnv 



4{6'u62)- 



Table [2] shows the notation we for each test statistic, along 
with the values of v = [h' {l)]'^ . 

We are now in position to state the following result. 



Proposition \_Let m and n be large, and consider the out- 
come 5'^(0i, 02) = s. The null hypothesis 0i = 02 can be 
rejected at a nominal level rf ifVr[xM ^ ^) —V- 

3.2 Kolmogorov-Smirnov contrast measure 

The empirical function of the sample Z = (Zi , . . . ,Z„) of a 
random variable Z is given by Fz{z) = n~^#{j : Zj < z}. 
Empirical functions can be used to test the null hypothesis 
Hq : Fzf = Fz" against the alternative H\ : Fz' ^ Fz" with ran- 
dom samples Z' and Z" of sizes m and n, respectively, using 
the test statistic ll44ll 



%s(z',z'0 = 



-— - max \Fz'{z)-Fz"{z)\, 



where Jks is the empirical estimate for the Kolmogorov- 
Smirnov distance Jks- In this case, the analytic expression 
for Jks under distributions with same number of looks is 
given by 



Jks(Z',Z'0= max 

— oo<Z<oo 

r(L-ai) 



L-ai,L+l, 



L-a2,L+l, 



-Lz 
71 

-Lz 

72 



where 2F1 is the hypergeometric function. 

Proposition 2 Assume Sks{^' ^^") = s. The null hypothesis 
Fz' = Fz" can be rejected at level r\ //'Pr[5'Ks(Z',Z'') > 5] = 
1 -L{s) < 7], where L{s) = 1 -2i:^^i(-l)^-i exp(-2y^V) 
and 5 is a critical value given in I144\I , provided that m and 
n are large. 

In terms of image analysis. Propositions [T] and [2] offer 
methods to statistically refute the hypothesis that two sam- 
ples obtained in different regions can be described by the 
same distribution. Such assessment is relevant, for instance, 
when choosing input training samples in supervised learning 
procedures, or when proposing similarity -based techniques 
such as segmentation and classification. 

Since the aforementioned results are asymptotic and 
valid only when the underlying hypothesis holds, the fol- 
lowing section shows a Monte Carlo approach for their as- 
sessment with finite samples and under contamination. Such 
assessment is needed since finite samples and outliers are 
unavoidable in practical situations. 

4 Validation 

This section presents a qualitative assessment of the test pro- 
posed for checking whether two samples are distributed ac- 
cording to the same law when finite size samples are used. 
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Table 2 {h, (^)-distances, statistics, h and (j) functions, and related quantity 

-distance Statistic h{y) v = [/z'(0)(^"(l)]-i 

Kullback-Leibler (Jkl) ^S'kl y/2 (x-l)logx 1 

Triangular (dj) Sj y,0<y<2 1 

Bhattacharyya (Jb) -\og{-y + 1),0 <y < 1 -v^+^ 4 

Arithmetic-geometric (Jag) '^'ag J ^log^ + (^) 4 



The proposed tests are submitted to uncontaminated data in 



Section [4^1] Subsequently, a contamination model based on 
innovative outliers, inspired by corner reflectors, is consid- 
ered in Section l4~2l The effect of contamination on the ML 
estimates of parameters is also assessed. 



4.1 Analysis with Simulated Data 

Speckled data simulated with the law were employed 
to assess the discussed parametric and nonparametric test 
statistics. The situations herein assessed were three levels 
of roughness (a G {—1.5,— 3,— 5}), three mean return va- 
lues (ji e {1,2,5}), and two values of the number of looks 
(L e {1,8}). The mean ji = -//(I + a) is directly related 
to the image brightness. We performed the tests at the 1% 
significance level using square windows of 7 x 7, 9 x 9, and 
11x11 pixels. Therefore, the considered sample sizes were 
A^ = m = ^G {49,81,121}. 

Empirical size and power of the proposed tests were 
sought as a means to measure and compare their perfor- 
mances. Monte Carlo experiments under five different sce- 
narios were used. Considering two image regions specified 
by parameter vectors (ai,/ii,L) and (a2,/i2,^), the follo- 
wing scenarios were set: (i) ai = a2, /ii =1^2, (ii) cci = oci, 
All ^ Ai2, (iii) cci > 0^2, Ml > M2, (iv) ai > a2, Mi < M2, 
and (v) ai 7^ a2, /ii = /i2- Whereas the first situation cor- 
responds to image regions of same roughness and bright- 
ness, scenario (ii) accounts for images of same roughness 
but with different brightness. Images of both distinct rough- 
ness and brightness are captured by scenarios (iii) and (iv). 
Additionally, the fifth scenario considers regions of different 
roughness and equal brightness. 

Fig. |4] shows the mean empirical test sizes, i.e., under 
scenario (i). With the exception of Sj, significance levels are 
above the nominal value (shown as a dashed line). The tests 
tend to reject more often when a increases, with little effect 
from sample sizes. In general, Sj is the most efficient mea- 
sure regarding the test size, since its estimative is in mean 
the closest to nominal values; Sj is immediately followed by 
^KS • As will be seen later, 5'ks surprisingly loses this ability 
in the presence of contamination. 

Goudail and Refregier | 22 | present graphically a rela- 
tionship between test power and the Kullback and Bhatta- 
charyya distances. Similarly, we illustrate the estimates of 



test powers with respect to the discrepancy measure given 
by 

5=(<iKL+^T+^B +t/AG+^Ks)/5, 

where, for the sake of simplicity, arguments (0i,02) were 
suppressed. 

Tables |3][6] present the test empirical power in scenar- 
ios (ii)-(v), respectively. As a general behavior, the empiri- 
cal power increases with the number of looks (L) and/or the 
sample size (N). 

Considering a fixed roughness. Table [3] provides evi- 
dence that the brightness change has more influence on the 
power in homogeneous regions. Fig. [5] illustrates the results 
of this table with respect to D for different sample sizes. In 
this case, we noticed that parametric tests outperformed the 
Kolmogorov-Smimov test. 

In Tables [4j [5] and |6j the inequality Power(5'Ks) < 
Power(^T) < Power(^B) < Power(^KL) < Power(^AG) is 
satisfied in every case of scenarios (iii) and (iv). For sce- 
nario (v), this fact holds when the number of looks is 8. 
Fig. [5] presents the powers as functions of D for different 
roughness pairs. The power increases with the sample size 
and with 5, as expected. These results show that the Sag test 
outperforms other techniques with respect to the test power, 
since it is able to discriminate more often samples from dif- 
ferent distributions. 



4.2 ML Estimation in Contaminated Data 

In practice, it is quite hard to guarantee that all observations 
in a sample come from random variables with the same dis- 
tribution. When collecting samples by visual inspection, fre- 
quently the onset of some type of contamination is experi- 
enced. 

One of the most common sources of contamination in 
SAR imagery is related to the double bounce phenomenon: 
a few pixels exhibit a very high return caused by either nat- 
ural or man-made targets that send back most of the energy 
they receive. This is the case of flooded forests and urban 
areas. The presence of a single pixel suffering from double 
bounce may produce meaningless estimates and test statis- 
tics, unless they exhibit some robustness with respect to this 
kind of contamination. Fig. [7] shows an example of the dou- 
ble bounce effect on an homogeneous field. 
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0,5 1.0 1.5 0.5 1.0 1.5 



Fig. 4 Mean of empirical test size scaled by 10 at 1% nominal level. 



Table 3 Rejection rates of {h, (^) -divergence and Kolmogorov-Smimov tests for scenarios (ii) 



M2/M2 


L N 


oc Sks 




St 


Sb 


^AG 


oc Sks Skl Sj Sb Sag 


« Sks Skl Sj Sb Sag 


2.0 


1 49 
81 
121 


-1.5 20.58 
48.36 
68.85 


27.84 
53.24 
75.04 


22.11 
48.79 
72.68 


26.89 
52.51 
74.84 


30.85 
55.52 
75.85 


-3.0 30.31 39.34 30.70 37.87 42.63 
61.85 68.03 63.52 67.26 70.14 
82.13 88.64 86.89 88.42 89.42 


-5.0 34.16 43.97 34.89 42.58 48.26 
68.64 75.65 70.50 75.01 77.89 
86.33 94.30 92.73 93.93 94.58 




8 49 
81 
121 


60.40 
91.13 
98.55 


78.53 
96.20 
99.73 


72.00 
94.65 
99.62 


78.02 
96.04 
99.73 


81.15 
96.69 
99.75 


90.85 97.41 95.89 97.34 97.78 
99.75 99.84 99.82 99.84 99.84 
100.00 100.00 100.00 100.00 100.00 


98.40 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


2.5 


1 49 
81 
121 


42.51 
77.22 
93.53 


53.07 
82.18 
95.73 


46.66 
79.73 
95.13 


52.46 
82.05 
95.68 


56.36 
83.41 
96.08 


56.80 70.15 62.23 69.00 73.33 
89.64 93.87 92.03 93.78 94.38 
98.27 99.54 99.40 99.50 99.56 


64.00 76.25 68.30 75.20 79.03 
93.47 97.02 95.72 96.86 97.42 
99.24 99.93 99.88 99.93 100.00 




8 49 
81 
121 


88.89 97.16 95.65 97.13 97.85 
99.53 99.95 99.91 99.95 99.96 
100.00 100.00 100.00 100.00 100.00 


99.60 99.96 99.95 99.96 99.98 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


99.96 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


5.0 


1 49 
81 
121 


96.58 98.91 98.27 98.95 99.06 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


99.29 98.55 97.39 98.42 98.81 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


99.76 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 




8 49 
81 
121 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 



In order to assess the robustness of the considered tests, 
we assume contamination in a parametric way: the data in 
a sample come either from the hypothesized distribution 
with probabiHty 1 — e, or from a scaled version of the same 
law with probabiHty £. The considered scaled distribution is 
such that its mean value is one hundred times larger than the 
mean of the original distribution. Let A be the event "pres- 
ence of outlier", which has Pr(A) = e, and A' its comple- 
ment. By the total probability rule, the cumulative distribu- 
tion function of the return is 



Following |[T]|5l, we chose to assess the situations without 
contamination e = and with mild levels of contamination 
{10-4,5x10-3}. 



Fz{z)=Pr{Z<z) 

= Pr ({Z < z} HA') + Pr ({Z < z} HA) 

= Pr(A')Pr(Z < z \ A') +Pr(A)Pr(Z <z\A) 



(l-^)^^0(a,r,L)W + 



(a,1007,L) 



Fig. 7 A comer reflector producing double bounce on a homogeneous 
background 
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Table 4 Rejection rates of {h, (^) -divergence and Kolmogorov-Smimov tests for scenarios (iii)-(v) for roughness pair (ai , = (—1.5, —3) 



Constants y Mi > M2 7 Mi = M2 7 Mi < M2 



L ( 71,72) ^KS ^KL ^AG ( 7l,72) ^KS ^KL ^AG {YU^ ^KS ^KL ^B ^AG 

1 49 (1,2) 1.05 7.12 1.17 3.53 11.94 (0.5,2) 13.45 16.15 11.64 15.14 18.56 (0.5,4) 82.27 89.51 86.89 89.48 90.47 

81 2.47 14.10 3.45 9.46 20.33 31.80 35.39 30.38 34.11 38.25 98.53 99.28 99.20 99.26 99.35 

121 2.64 24.98 10.20 19.17 31.83 51.13 58.36 54.75 57.68 60.30 99.89 100.00 100.00 100.00 100.00 



8 49 2.58 29.18 12.68 23.52 36.20 62.35 92.41 87.65 91.71 94.06 99.98 100.00 100.00 100.00 100.00 

81 7.33 50.77 36.49 47.39 56.12 93.31 99.51 99.16 99.45 99.53 100.00 100.00 100.00 100.00 100.00 

121 12.40 73.73 64.55 71.93 77.13 99.22 99.85 99.84 99.84 99.85 100.00 100.00 100.00 100.00 100.00 



1 49 (2.5,2) 69.45 89.53 82.24 87.95 91.89 (1,4) 13.78 15.55 10.85 14.41 18.28 (0.5,10) 99.95 100.00 100.00 100.00 100.00 

81 96.51 99.24 98.39 99.17 99.45 31.65 35.47 30.65 34.15 37.74 100.00 100.00 100.00 100.00 100.00 

121 99.62 100.00 99.98 100.00 100.00 49.98 58.51 54.67 57.57 60.48 100.00 100.00 100.00 100.00 100.00 



8 49 99.49 100.00 100.00 100.00 100.00 62.35 92.19 87.23 91.47 94.03 100.00 100.00 100.00 100.00 100.00 

81 100.00 100.00 100.00 100.00 100.00 92.93 99.55 99.29 99.49 99.62 100.00 100.00 100.00 100.00 100.00 

121 100.00 100.00 100.00 100.00 100.00 99.47 99.98 99.96 99.98 100.00 100.00 100.00 100.00 100.00 100.00 



1 49 (2.5,4) 5.29 18.77 6.70 13.06 24.88 (2.5,10) 13.71 14.83 10.65 14.00 17.47 (1,10) 94.65 97.74 96.96 97.74 97.92 

81 12.44 37.45 19.88 31.21 44.52 32.27 33.52 28.79 32.34 36.47 99.93 100.00 99.98 100.00 100.00 

121 20.35 57.55 41.69 53.04 63.38 50.87 58.23 54.08 57.37 60.09 99.98 100.00 100.00 100.00 100.00 



8 49 18.71 51.03 33.57 46.56 57.55 62.78 92.53 87.72 91.97 94.00 100.00 100.00 100.00 100.00 100.00 

81 44.84 78.12 67.69 76.02 81.37 92.11 99.58 99.25 99.56 99.69 100.00 100.00 100.00 100.00 100.00 

121 66.75 93.38 90.38 92.84 94.44 99.27 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 



Table 5 Rejection rates of {h, (^) -divergence and Kolmogorov-Smimov tests for scenarios (iii)-(v) for roughness pair (ai , = (— 1-5, —5) 



Constants 7 




Ml >M2 




7 




Ml = M2 7 


Ml < M2 


L 


(71,72) 






St 


Sb 


^AG (71,72) 


Sks 


Skl St Sb Sag (71,72) 


^KS ^KL ^B ^AG 


1 


49 (1,4) 
81 
121 


0.65 
1.33 
1.31 


11.15 
25.82 
44.16 


0.85 
3.85 
15.45 


4.07 
13.12 
32.35 


19.42 (0.5,4) 

35.38 

53.23 


25.36 
54.31 
75.85 


29.40 21.80 27.09 33.07 (0.5,8) 
60.22 53.02 57.79 63.30 
84.22 81.42 83.47 85.89 


91.91 95.52 94.08 95.55 95.96 
99.76 99.92 99.92 99.92 99.92 
100.00 100.00 100.00 100.00 100.00 


8 


49 
81 
121 


3.02 
11.13 
20.38 


74.90 
93.72 
97.87 


55.05 
88.57 
97.43 


69.6 
92.61 
97.76 


80.11 
94.8 
98.03 


92.13 99.91 99.77 99.89 99.96 
99.95 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


1 


49 (2.5,4) 
81 
121 


60.47 90.83 78.27 87.63 93.55 (1,8) 
91.87 99.48 98.51 99.32 99.66 
98.84 100.00 100.00 100.00 100.00 


24.11 
53.78 
75.71 


29.42 21.33 26.57 33.43 (0.5,20) 
59.83 53.71 58.03 62.63 
84.85 81.74 84.20 86.11 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


8 


49 
81 
121 


99.49 100.00 99.98 100.00 100.00 
99.98 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


92.13 99.89 99.79 99.89 99.91 
99.80 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


1 


49 (2.5,8) 
81 
121 


2.18 
7.02 
10.25 


21.98 
45.56 
68.37 


4.46 
15.82 
42.64 


11.99 
33.52 
59.95 


31.40 (2.5,20) 

55.15 

75.50 


23.85 
54.78 
75.38 


27.70 21.51 25.43 31.72 (1,20) 
60.56 53.55 58.22 63.61 
85.19 81.98 84.22 86.58 


97.49 98.86 98.67 98.83 99.02 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


8 


49 
81 
121 


13.67 
38.49 
62.85 


79.52 
96.08 
99.69 


58.63 
91.27 
99.05 


74.64 
95.11 
99.6 


84.73 
96.94 
99.84 


92.73 99.89 99.74 99.87 99.91 
99.89 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 



Simulation results are based on 5500 independent repli- 
cations, and instances were considered valid only if a G 
[10a, a/20], i.e., censoring was applied as in |15,40|. Ta- 
bles [7] [8j and |9] present the estimated bias (B), the coeffi- 
cient of variation (CV, in percentage), and the mean squared 
error (MSE) of the ML estimators a and 7 for the roughness 
degrees —1.5, —3, and —5 according to several simulation 
conditions and contamination levels, respectively. 

Let us first consider the bias for uncontaminated data. 
The only situations for which the bias reduces with increa- 
sing sample size are: (a) a = — 1.5, (b) a = — 3, L = 8, and 
(c) a = — 5, L = 8. This seemingly surprising situation is 
frequent when dealing with the distribution. Indeed, the 
information content is reduced in homogeneous and hetero- 
geneous situations; therefore, the good asymptotic proper- 



ties of ML estimators are not verified in samples of the con- 
sidered size |[6l [T5ll46l . An increase of the number of looks 
can partially compensate this behavior. 

This loss of information content is also verified by the 
dependence of the bias on the number of looks. It is ex- 
pected that the former is reduced increasing the latter, but 
this is only verified when a = —1.5, i.e., in extremely het- 
erogeneous areas. When dealing with heterogeneous data 
(a = —3), this reduction only takes place starting at A/^ = 81, 
and with homogeneous data (a = —5) at N = 121. Again, 
the more homogeneous the sample is, the smaller its infor- 
mation content becomes under the model. Regarding the 
N = 121 and L = 8 situation, which is expected to reveal 
the asymptotic behavior of the estimators, the bias increases 
when a is reduced from —1.5 to —5. 
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Table 6 Rejection rates of {h, (^) -divergence and Kolmogorov-Smimov tests for scenarios (iii)-(v) for (ai , = (—3, —5) 



Constants 7 




Ml >M2 






7 


Ml = M2 7 


Ml <M2 




^ (71572) 








^AG (71,72) 


Sks Skl Si 5b ^ag (71,72) 


^KS ^KL ^AG 


1 


49 (4,4) 
81 
101 

IZl 


19.2 
44.62 
66.18 


32.15 20.93 
62.18 51.78 
85.15 80.11 


28.45 
59.56 
84.24 


37.20 
65.80 
86.91 


(2,4) 


1 07 58 39 54 1 00 S^l 

L .\j 1 \j.uo yj.jyy yj.u^ y^L^o) 

1.96 0.98 0.58 0.82 1.49 
2.13 1.52 0.89 1.12 2.56 


46 4 S S9 SI 45 57 5180 S5 28 

'-ryj.'-ru uZ^.uL '-ru.ui uL.Oy uu.jLO 

81.96 83.39 80.52 83.33 84.91 
95.53 97.25 96.89 97.22 97.48 


8 


49 
81 
121 


81.25 
98.49 
99.98 


94.27 91.28 
98.75 98.6 
98.21 98.2 


93.94 
98.75 
98.21 


95.16 
98.79 
98.21 




^71 17 14 6 ID 70 IS 60 
9.22 26.17 19.21 24.48 29.42 
14.55 42.30 36.07 41.31 45.10 


QQ f^'y QQ Q7 QQ Of) QQ Q7 QQ Q6 
yy.Kjz^ yy.yz^ yy.y\j yy.yz^ yy.y\j 

100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


1 


49 (10,4) 
81 
121 


98.71 99.81 99.61 99.77 99.81 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


(4,8) 


1.16 0.54 0.38 0.50 1.15 (2,20) 
1.69 1.36 0.81 1.11 1.90 
2.53 2.11 1.27 1.71 2.93 


99.78 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


8 


49 
81 
121 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 




3.15 13.05 6.36 11.37 16.25 
8.56 25.55 18.48 23.98 28.87 
12.98 43.64 36.85 42.35 46.68 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


1 


49 (10,8) 
81 
121 


46.93 
81.65 
94.95 


65.05 52.45 
91.53 87.37 

99.06 98.55 


62.71 
90.84 
98.88 


69.91 (10,20) 

92.98 

99.36 


0.91 0.70 0.31 0.54 1.05 (4,20) 
2.18 1.30 0.96 1.06 2.11 
2.29 2.19 1.25 1.58 3.10 


74.84 81.73 76.48 81.42 84.18 
96.71 98.04 97.77 98.07 98.25 
99.71 99.95 99.92 99.95 99.95 


8 


49 
81 
121 


99.22 99.98 99.94 99.98 99.98 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 




2.96 12.71 6.67 11.31 16.28 
8.91 26.29 19.33 24.93 29.79 
13.58 42.24 36.07 41.11 44.92 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 




Regarding the influence of contamination on the bias, 
in all but a single case (namely A/^ = 49,L = l,e = 10~^), 
increasing the former increases the latter for both estimators. 
This increase in the bias ranges from, approximately, 34% 
(a = -1.5,L=l,7=0.5,£ = 10-'^)tolOOO%(a = -3,L = 
1 , 7 = 2, £ = 5 X 10"^), in both cases = 49. 

The coefficient of variation is reduced by increasing L 
and/or N and/or a, as expected, in uncontaminated data. In 
the presence of contamination, this behavior is lost with res- 
pect to a. 

The usual behavior of reduced MSE with increasing L 
and/or N and/or a is observed in the absence of contamina- 
tion. The presence of contamination consistently increases 
the MSE, except in four out of 324 situations, where ran- 
dom fluctuations masked the results. For instance, the MSE 
of 7 doubled in the a = -5, 7 = 20, L = 1, A/^ = 49 situation 
when e = 10~^, i.e., the mildest level of contamination. If 
a = -5, 7 = 20, L = 8, = 121, and e = 5 x 10"^ i.e., the 
strongest contamination, the MSE was affected by a factor 
of 13794. 



This strong influence of contamination on the MSE 
warns users about the need of quality samples and of robust 
estimators when dealing with this type of data. 

4.3 Test Sizes in Contaminated Data 

Based on the above results, the forthcoming robustness 
study is limited to the best statistics, namely 5't and 
along with which is expected to be robust. For these 
measures. Fig. [8] shows the behavior of the test size in the 
presence of fixed levels e of contamination. 

Under contamination, larger windows imply better em- 
pirical test sizes for the parametric tests. Tests based on 5't 
are the most robust for homogeneous regions. The empiri- 
cal test sizes of (/z, 0) -distance based tests for a = —1.5 are 
greater than the ones observed when aG{— 3,— 5}; this be- 
havior is more pronounced as e increases. 

Unexpectedly, the Kolmogorov-Smimov test tends to be 
less robust than parametric tests as contamination increases. 
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Fig. 6 Empirical test power in situations (iii)-(v) for different contrast measure at 1% nominal level. 



Fig. [8] indicates a low variability of the based test size 
across different simulation specifications. This suggests that 
image roughness and brightness do not significantly affect 
these particular estimated test sizes. 



5 Results 

In the following we present an application of the proposed 
tools to real data. The performance of the considered hy- 
pothesis tests are quantified by means of their empirical 
sizes and powers. Based on Fig.|9] samples of same and dif- 
ferent classes are constrasted. 

The image was obtained by the E-SAR sensor over sur- 
roundings of Munich, Germany [ 25 1, and its estimated num- 
ber of looks is 3.2. The area exhibits three distinct types of 
target roughness: (i) homogeneous (corresponding to pas- 
ture), (ii) heterogeneous (forest), and (iii) extremely hetero- 
geneous (urban areas). Samples were selected and submitted 
to statistical analysis, after parameter estimation by maxi- 
mum likelihood methods. 

Table 10 shows the a, 7, and /I estimates (whose in- 



sample size (# pixels). Values of /i for urban regions are in- 
finite in accordance to the theoretical background as seen in 
Equation (|4]). The table also presents the number of disjoint 
7x7 pixels blocks (# parts) which were used to estimate the 
test sizes and powers. 

The estimated roughness and brightness parameters sat- 
isfied the following inequalities: Ofpasture-i < ^asture-2 < 

^asture-3 ^ forest ^ ^rban-3 ^ ^rban-1 ^ ^rban-2 and 
Mpasture-1 ^ Mpasture-2 ^ Mpasture-3 ^ Mforest ^ Murban-3- Those 

inequalities confirm the interpretability of the parameters, 
and show how hard it is to discern different samples from 
the same class. 



Table 10 Parameter estimates 
Regions a y 



# pixels # parts 



terpretability was discussed in Section [471]), along with the 



pasture- 1 


-15.702 


39259 


2670.32 


1235 


25 


pasture-2 


-12.698 


80320 


6866.13 


1216 


24 


pasture-3 


-11.304 


162292 


15750.39 


1602 


32 


forest 


-9.339 


661183 


79288.04 


1606 


32 


urban- 1 


-0.759 


148413 


00 


2005 


40 


urban-2 


-0.388 


110183 


00 


3481 


71 


urban-3 


-1.079 


55583 


703582.30 


4657 


95 
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Table 7 Maximum likelihood estimates bias (B), coefficients of variation (CV%), and mean squared errors (MSE) for high roughness 



a Y L £ 



-1.5 0.5 1 
8 

1 1 
8 

2.5 1 

8 



N = 49 
B CV% MSE 



0.23 37.64 0.47 
0.11 23.58 0.15 
0.22 37.70 0.50 
0.10 22.73 0.15 
0.20 36.63 0.47 
0.12 23.31 0.15 



a 



N = U 
B CV% MSE 



0.18 32.50 0.33 

0.06 17.39 0.08 

0.17 32.24 0.34 

0.06 17.30 0.08 

0.19 32.59 0.33 

0.06 17.38 0.08 



0.27 32.33 0.40 

0.12 18.03 0.10 

0.26 31.36 0.37 

0.13 17.45 0.10 

0.28 32.05 0.40 

0.12 17.54 0.09 



0.84 42.28 1.68 

0.38 19.51 0.28 

0.82 40.98 1.58 

0.38 19.06 0.27 

0.85 41.88 1.69 

0.37 18.84 0.26 



1.42 40.45 3.42 
0.69 20.99 0.68 

1.43 39.83 3.40 
0.68 20.33 0.66 
1.42 40.18 3.40 
0.68 20.17 0.65 



A^= 121 
B CV% MSE 



0.14 27.15 0.22 

0.04 13.95 0.05 

0.12 26.86 0.21 

0.05 14.01 0.05 

0.13 26.74 0.20 

0.04 13.87 0.05 



0.23 27.15 0.27 
0.10 14.08 0.06 
0.23 27.19 0.28 
0.10 14.16 0.06 
0.23 26.78 0.26 
0.10 14.20 0.06 



0.75 35.29 1.19 

0.35 15.16 0.20 

0.75 34.74 1.17 

0.35 15.14 0.20 

0.75 34.69 1.16 

0.35 15.12 0.20 



1.41 36.47 3.11 
0.64 16.35 0.53 

1.42 36.66 3.16 
0.64 16.13 0.53 
1.41 36.94 3.16 
0.64 16.01 0.53 



A^ = 49 
B CV% MSE 



0.12 56.10 0.13 

0.05 29.83 0.03 

0.24 56.60 0.60 

0.09 28.24 0.11 

0.55 54.94 3.47 

0.26 29.26 0.66 



0.19 52.79 0.17 

0.09 28.50 0.04 

0.37 53.47 0.67 

0.18 28.72 0.15 

0.91 52.60 4.05 

0.47 28.79 0.95 



0.74 60.30 1.10 

0.35 29.15 0.19 

1.44 60.88 4.30 

0.70 29.95 0.75 

3.65 60.33 27.04 

1.76 29.81 4.73 



1.93 53.58 5.41 

1.05 29.32 1.30 

3.84 55.04 21.86 

2.10 30.80 5.31 

9.69 55.03 138.98 

5.28 30.76 33.60 
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A^ = 81 
B CV% MSE 



0.09 47.64 0.09 

0.03 21.76 0.01 

0.19 47.66 0.37 

0.05 21.78 0.06 

0.49 48.49 2.25 

0.14 21.88 0.35 



0.16 46.73 0.12 

0.07 21.98 0.02 

0.31 45.41 0.45 

0.14 20.98 0.08 

0.83 45.77 3.01 

0.34 21.54 0.49 



0.66 53.66 0.82 

0.32 22.41 0.14 

1.31 52.46 3.17 

0.64 21.67 0.53 

3.34 52.76 20.68 

1.58 21.76 3.28 



1.87 48.52 4.82 

1.00 23.29 1.11 

3.76 47.85 19.30 

1.98 22.65 4.37 

9.36 48.42 120.59 

4.95 22.63 27.39 



A^= 121 
B CV% MSE 



0.07 40.00 0.06 

0.02 17.67 0.01 

0.13 39.59 0.23 

0.04 17.64 0.04 

0.35 39.55 1.35 

0.09 17.45 0.22 



0.14 39.15 0.08 

0.06 17.14 0.01 

0.28 38.05 0.31 

0.12 17.53 0.05 

0.67 38.11 1.91 

0.31 17.48 0.33 



0.60 44.81 0.60 

0.30 17.38 0.11 

1.20 43.15 2.33 

0.61 17.52 0.45 

2.98 43.82 14.68 

1.52 17.47 2.81 



1.85 43.66 4.45 

0.96 18.29 1.00 

3.71 43.21 17.90 

1.92 18.05 3.96 

9.24 43.82 111.92 

4.82 17.80 24.89 



10 



10 



5x 10 



0.31 36.75 0.54 

0.17 23.78 0.19 

0.29 36.41 0.51 

0.17 23.51 0.18 

0.30 36.71 0.52 

0.17 23.34 0.18 



0.94 47.42 2.22 

0.44 25.90 0.44 

0.91 47.46 2.14 

0.44 25.98 0.44 

0.93 48.01 2.23 

0.44 25.83 0.45 



1.45 44.24 3.79 
0.75 26.83 0.93 
1.43 44.82 3.78 
0.75 27.64 0.95 

1.46 45.50 3.93 
0.76 27.45 0.97 



Table 8 Maximum likelihood estimates bias (B), coefficients of variation (CV%), and mean squared errors (MSE) for mid roughness 



a Y L £ 



-3.0 2.0 1 
8 

4.0 1 

8 

10.0 1 
8 



A^ = 49 
B CV% MSE 



0.15 42.11 1.81 

0.31 27.28 0.90 

0.17 41.35 1.79 

0.29 27.51 0.88 

0.18 41.92 1.84 

0.28 27.79 0.88 



a 



N = Sl 
B CV% MSE 



0.29 37.21 1.62 
0.18 22.14 0.49 
0.29 37.47 1.63 
0.18 21.29 0.49 
0.32 38.17 1.66 
0.16 21.15 0.51 



0.43 37.48 1.84 
0.29 21.30 0.57 
0.47 37.01 1.87 
0.29 21.72 0.59 
0.40 36.90 1.74 
0.30 21.61 0.60 



1.14 31.11 2.95 

1.03 21.99 1.84 

1.16 30.90 3.01 

1.03 22.05 1.84 

1.10 31.03 2.84 

1.05 22.39 1.93 



1.93 24.11 5.12 

2.18 18.18 5.62 

1.96 23.49 5.18 

2.18 17.96 5.64 
1.88 23.62 4.86 

2.19 18.12 5.69 



A^= 121 
B CV% MSE 



0.29 34.61 1.42 
0.12 17.58 0.30 
0.30 34.99 1.45 
0.11 18.28 0.32 
0.31 35.07 1.51 
0.12 17.93 0.29 



0.46 34.04 1.60 
0.23 17.48 0.37 
0.47 33.76 1.59 
0.23 17.18 0.36 
0.48 33.94 1.63 
0.23 17.03 0.36 



1.30 28.34 3.18 
0.97 18.74 1.50 

1.34 27.64 3.22 
0.97 18.47 1.48 

1.35 28.24 3.33 
0.99 18.49 1.52 



2.19 20.16 5.91 
2.28 15.81 5.87 

2.26 19.92 6.20 

2.27 15.97 5.88 
2.14 21.23 5.76 
2.30 15.95 6.00 



A^ = 49 
B CV% MSE 



0.16 56.46 1.55 

0.24 31.38 0.55 

0.35 55.24 5.95 

0.45 31.65 2.14 

0.90 55.33 39.62 

1.09 31.91 13.19 



0.33 54.06 1.69 

0.35 31.15 0.66 

0.75 52.75 6.85 

0.68 30.78 2.53 

1.68 54.01 42.65 

1.68 30.98 15.93 



1.22 44.44 3.55 

1.27 28.36 2.49 

2.70 43.67 15.85 

2.55 28.48 10.00 

6.31 44.07 91.51 

6.35 28.37 61.86 



3.99 33.89 20.00 
4.05 23.00 18.33 
7.85 34.20 77.99 
8.10 22.78 73.24 
19.13 35.19 471.04 
20.35 22.38 460.41 
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A^ = 81 
B CV% MSE 



0.25 48.36 1.27 

0.14 25.21 0.29 

0.53 48.57 5.35 

0.28 24.47 1.17 

1.42 49.42 32.20 

0.65 24.30 7.71 



0.42 48.17 1.54 

0.26 24.49 0.37 

0.92 47.45 6.29 

0.53 24.77 1.53 

2.05 47.10 36.42 

1.39 24.67 9.82 



1.49 38.04 3.98 
1.24 24.24 2.15 
3.04 37.61 16.23 
2.47 24.07 8.53 
7.47 38.16 100.21 
6.31 24.51 55.79 



4.33 28.75 22.08 
4.22 19.88 19.32 
8.79 28.54 90.59 
8.47 19.56 77.78 
21.52 27.91 540.58 
21.22 19.70 488.21 



A^= 121 
B CV% MSE 



0.25 44.09 1.12 

0.09 20.19 0.18 

0.53 44.84 4.33 

0.18 20.84 0.75 

1.38 44.93 29.69 

0.48 20.61 4.36 



0.44 43.13 1.31 

0.22 19.85 0.24 

0.91 43.08 5.30 

0.43 19.70 0.94 

2.31 42.67 32.91 

1.09 19.26 5.74 



1.65 34.06 4.28 
1.18 20.27 1.82 
3.35 33.24 17.21 
2.37 20.11 7.26 
8.47 33.70 110.56 
5.99 19.97 46.09 



4.72 23.54 24.78 
4.34 17.00 19.96 
9.61 23.78 102.92 
8.67 17.21 79.98 
23.29 25.41 614.05 
21.85 17.19 507.52 



10 



10 



0.32 41.92 2.03 
0.40 27.47 1.03 
0.33 40.13 1.90 
0.38 27.10 0.99 
0.3141.12 1.95 
0.39 27.33 1.01 



0.84 35.35 2.54 
1.07 25.91 2.24 
0.93 34.41 2.68 
1.07 25.88 2.25 
0.87 35.12 2.59 
1.07 25.95 2.26 



5 X 10^ 



1.64 27.29 4.29 
2.03 20.87 5.22 
1.62 28.30 4.34 
2.03 20.69 5.23 
1.52 28.95 4.02 
2.05 20.42 5.27 



Table 1 1 presents the rejection rates at significance levels 
1% and 10% of samples from the same area, i.e., the test size 
or Type I error. Such nominal levels were considered useful 
in practical applications (cf. reference (43J.) These results 
show that ^kl, and Sj tests have excellent performance 
with respect to this criterion, and that the hardest situation is 
related to urban areas. Furthermore, the 5'ag test shows op- 
timal size test for 1% nominal level, failing on urban- 1 and 



urban-2 situations at 10%, while the Sb test size is inappro- 
priate in case urban-2. 



Table [12] shows the test powers at 1% and 10% nom- 
inal levels, i.e., the inability to reject samples from differ- 
ent types. The best performances are observed in 
divergence tests, mainly when confronting areas whose dis- 
tributions differ more significantly; among them, Sj and 5'kl 
are the weakest ones. The worst result was observed in the 
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Table 9 Maximum likelihood estimates bias (B), coefficients of variation (CV%), and mean squared errors (MSE) for low roughness 



N = A9 
B CV% MSE 



A^ = 81 
CV% MSE 



B 



A^= 121 
CV% MSE 



0.22 39.77 
0.13 26.71 
0.32 40.23 
0.11 27.77 
0.18 39.53 
0.13 27.63 



0.39 37.67 
0.46 21.34 
0.49 38.28 
0.44 20.82 
0.48 38.78 
0.45 20.83 



4.27 
1.28 
4.23 
1.28 
4.13 
1.30 



4.26 
1.57 
4.66 
1.48 
4.70 
1.49 



1.73 30.71 
1.93 20.90 
1.76 30.81 

1.92 20.74 
1.71 30.66 

1.93 20.86 



7.29 
5.81 
7.45 
5.76 
7.16 
5.82 



B 



A^ = 49 
CV% MSE 



A^ = 81 
CV% MSE 



A^= 121 
CV% MSE 



0.30 58.61 
0.47 33.62 
0.48 57.82 
0.92 34.18 
1.09 58.27 
2.28 33.64 



4.73 
2.42 
19.89 
10.23 
122.69 
60.17 



0.02 57.34 
0.68 33.27 
0.05 56.87 
1.33 33.63 
0.37 57.22 
3.16 32.88 



5.21 
2.90 
20.43 
11.63 
126.31 
67.97 



1.48 46.67 
2.40 28.81 
3.05 46.15 

4.68 28.27 

7.69 47.46 
11.90 28.62 



8.71 
9.15 
35.33 
34.79 
231.81 
224.98 



a Y L £ 



-5.0 4.0 1 
8 

8.0 1 

8 

20.0 1 
8 



0.33 46.97 4.85 

0.52 30.38 3.01 

0.30 46.86 5.11 

0.51 30.85 3.14 

0.25 47.26 4.99 

0.52 30.61 2.99 



0.02 42.92 4.75 

0.35 28.62 2.04 

0.06 42.50 4.75 

0.28 29.04 2.08 

0.00 43.57 4.74 

0.31 28.25 2.04 



0.20 42.56 4.93 

0.57 25.28 2.31 

0.2141.15 4.65 

0.58 25.85 2.41 

0.17 41.88 4.71 

0.55 25.60 2.32 



1.38 34.95 6.87 

1.93 23.19 6.30 

1.38 33.55 6.48 

1.89 23.43 6.20 

1.32 33.96 6.36 

1.88 23.52 6.13 



2.66 26.17 29.23 
3.97 14.84 42.14 
2.74 24.23 32.54 
3.89 15.32 40.99 
2.64 23.50 25.52 
3.88 15.45 41.52 



3.39 19.56 41.10 
4.21 12.92 44.56 
3.19 22.56 43.42 
4.24 12.72 44.41 
3.26 21.64 36.30 
4.19 13.21 44.84 



5.69 35.86 100.86 
7.56 20.03 124.35 
12.41 33.54 441.66 
14.95 19.77 490.83 
30.16 34.71 1959.53 
37.37 20.75 3047.61 



0.03 51.94 4.58 

0.33 31.08 1.61 

0.05 51.74 17.98 

0.52 31.38 6.64 

0.05 53.04 110.26 

1.43 30.69 39.68 



0.27 51.54 4.93 

0.57 27.83 1.95 

0.58 49.38 18.29 

1.15 28.59 8.19 

1.15 50.32 114.65 

2.79 28.23 49.17 



2.10 40.65 10.57 

2.47 24.71 8.64 

4.26 39.08 41.11 

4.88 25.32 34.44 

10.21 39.33 245.49 

12.11 25.08 211.47 



6.88 30.02 115.62 
8.19 16.09 133.49 
13.86 27.47 498.89 
16.12 16.65 522.42 
34.69 27.82 2571.16 
40.32 16.61 3304.94 



0.22 47.82 3.90 

0.13 28.51 1.02 

0.64 48.68 15.97 

0.21 29.41 4.06 

0.89 47.03 98.05 

0.63 29.41 25.36 



0.45 44.61 4.15 

0.47 23.43 1.32 

1.09 45.65 18.43 

0.90 22.75 4.91 

2.62 45.63 113.44 

2.30 22.88 31.31 



2.49 35.23 11.44 

2.46 22.29 8.13 

4.94 35.14 45.13 

4.91 21.99 32.12 

12.31 35.57 283.43 

12.35 22.24 204.18 



7.95 22.52 151.55 
8.49 14.03 140.03 
15.01 25.64 632.03 
17.09 13.67 558.04 
38.22 23.54 3357.49 
42.35 14.20 3523.39 



10 



10 



0.14 46.39 
0.70 30.21 
0.14 46.19 
0.68 30.30 
0.19 46.55 
0.64 29.89 



5 X 10 



5.11 
3.44 
5.06 
3.43 
5.04 
3.26 



0.76 38.72 5.54 

1.84 26.58 6.68 

0.83 38.48 5.72 

1.80 26.35 6.47 

0.77 38.63 5.57 

1.84 26.76 6.72 



1.90 31.41 25.32 
3.52 18.38 38.51 
2.23 29.07 28.14 
3.46 18.40 38.25 
2.07 29.49 19.21 
3.45 19.05 37.62 




Table 11 Rejection rates of 
Hi: (ai,7i) ^ (a2,r2), with jii 
levels 



-divergence 
= jU2 at 1% and 



tests under 
10% nominal 



Fig. 9 E-SAR image and selected regions. 



urban- 1 vs. urban-2 situation whose samples are, in fact, 
very close and hard to discern. The Kolmogorov-Smirnov 
test produces unacceptable results when contrasting urban-3 
vs. forest, and urban- 1 vs. urban-2 samples. In particular, the 
situation urban-3 vs. forest contrasts two noticeably distinct 
regions, as shown in Fig.|9]and Table 10 Fig.[TO|presents the 
empirical distribution functions of the areas considered. It il- 





1% nominal level 




10% nominal level 




Regions 


^KS ^KL 


^AG 


^KS 




5't Sb 


^AG 


pasture- 1 
pasture-2 
pasture- 3 


0.33 0.00 0.00 
0.36 0.00 0.00 
0.60 0.00 0.00 


0.00 0.00 
0.00 0.00 
0.36 0.36 


5.00 
6.16 
6.05 


0.00 
0.00 
3.62 


0.00 0.00 
0.00 0.00 
3.62 3.62 


0.00 
0.00 
3.99 


forest 


2.22 0.00 0.00 


0.00 0.00 


9.68 


2.67 


1.67 2.00 


2.33 


urban- 1 
urban-2 
urban-3 


0.00 0.00 0.00 0.13 0.00 
0.00 0.00 0.00 99.96 0.00 
0.29 4.23 1.97 3.56 5.64 


2.18 
0.97 
5.98 


4.10 
2.01 
20.58 


2.18 34.49 90.64 
1.21 99.96 100.00 
16.84 20.85 31.02 



lustrates the close connection of the efficiency to these 
functions, regardless the texture differences, which makes 
the test an inadequate tool for discrimination in speck- 
led imagery. 



6 Conclusions 

This paper presented a comparison among the performance 
of four parametric tests based on stochastic distances and the 
Kolmogorov-Smirnov test. We also presented compact for- 
mulas for the Kolmogorov-Smirnov contrast measure, and a 
study of the robustness of the maximum likelihood estima- 
tor under the law. The assessment was made with Monte 
Carlo experiments varying the parameters of the distri- 
bution, regarded as a universal model for speckled data, and 
the level of innovative contamination. In doing so, we pro- 
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Fig. 8 Test sizes at 1% nominal level and £. For each sample size, all considered scenarios were lexicographically ordered: (/;a,/i) 
{(1; — 1.5, 1), . . . , (4; —3, 1), . . . , (7; —5, 1), . . .}]. Index / was employed as abscissa value. 
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Fig. 10 Empirical cumulative distribution function for different re- 
gions in SAR data for sample values 



vided new results regarding the both the accuracy and the 
robustness of the measures under examination. 



We show numerical evidence that the test based on the 
triangular distance has, in general, smaller empirical Ty- 
pe I error than the test based on the Kolmogorov-Smirnov 
distance 5'ks • Using a parametric and plausible contamina- 
tion model, we illustrate that when the chances of observing 
aberrant data increase, the test based on has its perfor- 
mance diminished with respect to the results obtained from 
the -distance methods. This is an important warning 
due to the widespread use of the Kolmogorov-Smirnov test 
when dealing with non-normal data. 

The 5'ag nieasure presented the best performance with 
respect to the test power. However, for a given number of 
looks, we observed that the test power performance of the 
four proposed parametric measures was roughly the same, 
and that they were more accurate than the test related to the 

nieasure. 

These conclusions suggest the use of the for SAR 
data processing due to its resulting accurate test empirical 
size. On the other hand, the S^g test is a more appropriate 
tool whenever the focus is on the test power. 
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Table 12 Rejection rates of {h, (^) -divergence tests under Hi : 7i) 7^ («2, 72), with /ii = /i2 at 1% and 10% nominal levels 





1 % nominal level 


10% nominal level 


Regions 


^KS ^KL ^AG 


^KS ^KL ^B ^AG 


pasture- 1 xpasture-2 
pasture- 1 xpasture-3 
pasture-2 xpasture-3 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 94.64 87.50 90.48 95.24 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 98.21 97.02 98.21 98.21 


forest X pasture- 1 
forest xpasture-2 
forest xpasture-3 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 


100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 



urban- 1 x pasture- 1 
urban- 1 xpasture-2 
urban- 1 xpasture-3 
urban-2 x pasture- 1 
urban-2 xpasture-2 
urban-2 xpasture-3 
urban-3 x pasture- 1 
urban-3 xpasture-2 
urban-3 xpasture-3 
urban- 1 X forest 
urban-2 X forest 
urban-3 X forest 
urban- 1 X urban-2 
urban-2 X urban-3 
urban- Ix urban-3 
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.40 
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.00 
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99 
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.17 
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We illustrated the discriminatory capability of these 
measures using real data. Tests based on the classical mea- 
sures here tested, namely the Kullback-Leibler, Kolmo- 
gorov-Smirnov, and Bhattacharyya distances, were outper- 
formed by the other distances and should not be employed. 

A Computational information 

Simulations of the (/z, (^) -distance based parametric tests were per- 
formed in Ox programming environment 1 10|; function QAGI was em- 
ployed for numerical integration. The Kolmogorov-Smimov test was 
coded in R and function ks.test was called |35|. The computation 
time of the triangular and arithmetic-geometric distances took typica- 
lly less than one and four millisecond, respectively, when performed 
in a Pentium processor at 3.20 GHz. We used the George Marsaglia's 
multiply-with-carry with 52 bits pseudorandom number generator, 
which has an approximate period of 2^^^^. 
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